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Abstract 

Large deviations is a framework for rigorously studying the probabilit- 
ies of rare events in terms of a rate function. We study the atom maser model 
in the context of large deviations, following the work done by Garrahan, Ar- 
mour and Lesanovsky [1 1 which relates this rate function to the occurrence 
of quantum dynamical phase transitions. We show that a large deviations 
principle holds for the associated classical counting process. This is done 
by deforming the Lindblad generator for the cavity dynamics and analys- 
ing the spectral properties of the associated deformed quantum dynamical 
semigroup. In particular, the large deviations rate function is given by the 
spectral bound of the generator. We use the rate function to discuss meta- 
stable behaviour of the stationary distribution of the counting process and 
its connection to quantum dynamical phase transitions. 

1 Introduction 

In this paper we study the atom counting statistics associated to the atom maser 
model. In particular, we are interested in whether any quantum dynamical 
phase transitions occur, using a large deviations approach. Previous work was 
done in P|j6j; in this paper, we expand on the mathematical background to 
the results obtained in The large deviations approach originated in non- 
equilibrium statistical mechanics |5|, where transitions between 'active' and 
'passive' dynamical phases were analysed using a large deviations rate func- 
tion. This rate function plays the role of a free energy in statistical mechanics, 
and this analogy is used to find dynamical phase transitions as non-analyticities 
in the rate function. 

In general, the large deviations approach involves quantifying rare events; 
in this case, the events are the trajectories followed by the state of the atom 
maser cavity. These trajectories can be characterised in term of the number of 



jumps between energy levels per unit time, and in this sense trajectories can be 
called active or passive. Trajectories which are highly active or highly passive 
are considered rare, and the large deviations approach is used to quantify their 
occurrence. The phase transitions we consider are between values of the dy- 
namical parameters for which active (or passive) trajectories suddenly become 
common. These points are closely related to the large deviations rate function, 
which quantifies the occurrence of such trajectories. 

In the atom maser model, atoms are prepared in a known state, and interact 
with a resonant cavity (which is in contact with a thermal bath) one at a time. 
We link properties of the stationary state of the cavity to the counting statistics 
of the output atoms. Our main result is that a large deviations principle holds 
for the counting process associated to atoms emerging from the cavity in the 
ground state. This means that quantum dynamical phase transitions do not oc- 
cur in the sense of non-analytic behaviour of the large deviations rate function. 
The proof of this is based on analysis of the spectral properties of the generator 
of the dynamics; in particular, we show that a Krein-Rutman theorem applies 
to the operators in the associated quantum dynamical semigroup. 

We note some other previous work on large deviations in quantum systems. 
In (7J (see also 18]) a large deviation principle is shown to hold for correlated 
states on quantum spin chains, a related (but finite-dimensional) system; large 
deviations for quantum Markov semigroups are studied in 1 9 1 . 

This paper is organised as follows. In Section 2 we introduce the back- 
ground to our problem, starting with the atom maser, where we discuss the 
master equation and the counting process. The rest of the section consists of 
the mathematical background, with one part about the theory of large devi- 
ations and the other part introducing some operator analysis. 

In Section 3 we state our main result, the counting process associated to the 
atom maser satisfies a large deviations principle; an outline of the proof is given. 
The proof of our result is found in Section 4. To show that the large deviations 
principle holds, we prove that a Krein-Rutman theorem holds for the generator 
of the atom maser dynamics. We discuss the results in Section 5, where we 
perform a numerical analysis of the problem. 

2 Background 

Atom maser 

The atom maser consists of a sequence of two-level atoms passing successively 
through a cavity and interacting resonantly with the electromagnetic field in- 
side the cavity. The two-level atoms enter the cavity in a fixed state, chosen here 
to be the excited state, and for simplicity we assume that only a single atom 
passes through the cavity at a time. In addition, the cavity is also coupled to 
a thermal bath which represents the interaction between the (non-ideal) cavity 
and the environment. The combined effects of the interactions with the atoms 
and the environment changes the state of the cavity whose time evolution is 
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described by a quantum Markov semigroup, in a certain coarse grained approx- 
imation; see Refs. |10| and [11 1 for a mathematical overview, and [12] for the 
physical derivation of the master equation. Here we will give an intuitive de- 
scription of the dynamics starting with a simplified discrete time model, with 
an emphasis on the statistics of measurements performed on the atoms. 

The cavity is described by a one mode continuos variables system with Hil- 
bert space t) = l 2 (N) whose canonical basis vectors (e n ) n ^ represent pure 
states of fixed number of photons. Therefore, if i|> G f) is a pure state, the photon 
number distribution of the cavity is given by |(e n ,i|j)j 2 . In general, (mixed) states 
are described by density operators i.e. positive trace-class operators p € T(f)) 
which are normalised to have unit trace, and the observables are represented 
by self -adjoint elements of the algebra of bounded operators S (F) ). Recall that 
the annihilation operator a on t) is defined by 



V^fxen-i ifn>0, 
ifn = 0, 



Its adjoint is the creation operator a* and N = a* a is the photon number operator. 
The atom is modelled on the two-dimensional Hilbert space C 2 with stand- 
ard orthonormal basis {|0), |1)} consisting of the 'ground' and 'excited' states. 
We denote by cr* and cr the corresponding raising and lowering operators (i.e. 
cr* |0) = 1 1 ) etc. ) The interaction between an atom and the cavity is described by 
the Jaynes-Cummings hamiltonian 

H int = ~9{ aa * + ff * a )- 

where g is the coupling constant which for simplicity is considered to be con- 
stant across the cavity. The free hamiltonian is 

Hjr ree = cua*a + a>cr*cr, 

where cu is the frequency of the resonant mode; however by passing to the in- 
teraction picture the effect of the free evolution can be ignored. Therefore if 
the interaction lasts for a time t , the joint evolution is described by the unitary 
operator U := exp(itoHj n |.) whose action on a product initial state is 

U : |k) ® |1) h> cos(4)Vk + l)|k) ® |1) + sin(cbVk + I)|k + 1) ® |0), 

where 4> := t g is the accumulated Rabi frequency. Conditional on the outcome 
of a measurement in the standard basis of the atom, the cavity remains in state 
|k) with probability cos 2 (4)v / k~+T) or it gains an excitation with probability 
sin (Wk+T). If we average over the outcome, we obtain a transfer operator 
T, acting on the cavity state space 

o;(p)=LipL; + L a pLS (i) 

where the operators Li and L2 are given by 

I_i = cos(4>Vaa*J, L 2 = a 7 = — . (2) 
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Since each atom interacts with the cavity only once, the state of the cavity after 
n such interactions is given by p(n) = 7^[p) which can be interpreted as a 
discrete time quantum Markov dynamics. 

Let A n denote the random variable equal to the number of atoms detected 
in the ground state after n interactions. The moment generating function (MGF) 
of A n is 

E (e s A " ) = Y_ p ( A n = k) e sk (3) 

where the expectation is taken with respect to the probability distribution of 
the counting process. Using <jT|) it is easy to see that the MGF is equal to 

E (e«An) =Tr(0? s (p))=Tr(p07(l)) 

where p is the inital state of the cavity and s (and its dual T s ) is obtained by 
biasing the second Kraus operator in Eq. Q by a factor e s : 

T* s (p) =LipLJ +e s L 2 pL*. 

As we will see, this allows us to study the probabilistic properties of the count- 
ing process A(t) in terms of spectral properties of the operator T s . Before that, 
we will make the model more realistic by passing to a continuous time descrip- 
tion in which the atoms are Poisson distributed in time with intensity N ex , and 
the cavity is in contact with a thermal bath. If one ignores the details of short 
term cavity evolution, the discrete time dynamics can be replaced by coarse 
grained continuous time Lindblad (master) equation 

ip(t)=£*(p(t)), 

4 



£*( P ) = 2l(L i pLt-i{LtL i) p}). (4) 



with jump operators Li defined by 




1-1 = V Nex COs(4)\/ QQ*), (5) 

, sin(4>V cm*) 



(6) 



L 3 = VV+la, (7) 
L 4 = ^a*. (8) 

As before, the operators and L 2 are associated to the detection of an atom 
in the excited, and respectively ground state. The emission and absorption of 
photons due to the contact with the bath is represented by the operators L3 
and respectively L4. The Heisenberg picture Lindbladian £ is the generator 
of a strongly continuous semigroup (cf . |13) for an introduction to the general 
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theory) on 23 (f) ). This means that there exists a family (T(t) ) t >o of maps on 23 ( f) ) 
which satisfy the semigroup property 



T(t)T(s) =T(t + s) forallt,s>0, 

no) = I, 

such that 1 1 — ^ T(t) (X) is norm continuous for all X e 23 (f) ) . The generator £ can 
be recovered by 

£(X)=limi(T(H)(X)-X), 

for all X in the domain of L . Although no simple expression exists for the op- 
erators T(t) in terms of the generator £, it can be helpful to think of T(t) as the 
exponential of the generator 

T(t)(X) = e t£ (X), (9) 

especially from the point of view of relating spectral properties of L to those 
of T(t), e.g. spectral mapping theorems. By definition, equation |9) does hold for 
analytic vectors of L which form a core of its domain. 

In the discrete-time model we were able to express the moment generating 
function in terms of a biased transition operator. In the continuous-time case, a 
similar expression can be obtained in terms of a biased generator £ s , where the 
jump operator L 2 associated to the ground state counting process is multiplied 
by the factor e s 

£,(X) = e s v 1 XL 1 -~{V 1 u,x}+'£_ (uxu-^atu,*}) • (io) 

Note that L s is still the generator of a completely positive semigroup, but it 
is no longer identity preserving, and therefore does not represent a physical 
evolution except for s = 0. Its importance lies in the fact that the moment gen- 
erating function associated to A t can be expressed in analogy to equation |3]l in 
the discrete time setting, as 

E(e sAt ) =Tr(pe t£s (l)) . (11) 

This equation plays a central role in this paper; we will use it to formulate a large 
deviations principle for the counting process A t which analyse the asymptotic 
(long-time) behaviour of the process. In particular, we will use Eq. (IT) to relate 
the moment generating function of A t to the spectral properties of £ s . 

In asymptotics it is important to verify whether the Markov semigroup is 
ergodic, i.e. there exists a unique stationary state p ss (i.e. £*(p ss ) = 0), such 
that 

lim T*(t)(p) = p ss , 

t— >oo 

for all initial states p. In the case of the atom maser, such a stationary state exists 
and is unique [14]. Since the semigroup generated by L leaves invariant the 
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commutative subalgebra generated by the number operator N, the stationary 
state must be diagonal in the number basis. The restriction of T(t) diagonal 
algebra we obtain a classical process known as a birth-death process, for which 
the (unique) stationary state can be computed explicitly as 

, , m n ( v N ex sin 2 (ctn/k]\ 
p ss (n) :== p ss (0) 11(^1+ — v ] • (12) 



k=l 



where p ss (n) is the probability of finding n photons in the cavity and p S s(0) is 
taken such that Tr(p ss ) = 1. The mean photon number of the stationary state is 
given by 

(N) = ^np ss (n). (13) 



The dependence of the mean photon number and photon number distribution 
on the parameter cf> is shown in Fig. [lj with the values v = 0.15 and N ex = 100 
for the other parameters. We note two interesting features in this figure: first, 
there is a sharp change in the mean photon number at cc = 1 followed by less 
pronounced jumps near a = 2n and a = 4tc. 

The other, related, feature to note concerns the photon number distribution, 
which only has a single peak for most values of a. In certain critical regions for 
a, such as 4.6 < a < 7.8, this distribution becomes bimodal, indicating that the 
stationary state is bistable. Such metastable behaviour of the atom maser has been 
investigated, using different techniques, in Ref. |15|. We study this metastable 
behaviour using the large deviations formalism. The large deviations rate func- 
tion, introduced below, can be used as an indicator of the sharp jumps in the 
mean photon number. 



Large deviations 

The main result of this paper is the existence of what is called a large deviations 
principle for the counting process A t introduced above. In 1 1 . 2] the stationary 
dynamics of the atom maser is analysed using the theories of large deviations 
and dynamical phase transitions. The large deviations approach raises import- 
ant questions related to the existence of dynamical phase transitions which can 
be formulated in terms of the spectral properties of the perturbed generator L s 
defined above. We note that a large deviations principle has been shown to hold 
in (7j in the context of quantum spin chains; the main difference is that we are 
looking at a continuous time Markov process on an infinite-dimensional space. 

We follow |lj|2j in using large deviations to analyse quantum dynamical 
phase transitions, where a biased generator is interpreted as creating traject- 
ories with a greater or lower number of ground state atoms emerging from 
the atom maser. This biased generator appears in Eq. [10] where we use it to 
define the moment generating function associated to the counting process A t . 
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Figure 1: Mean photon number (white line) and photon number distribution 
(background) in the stationary state p ss as function of a = -^/Nex^ 

This biased generator £ s , although it is not physical due to not being trace pre- 
serving, changes the expected number of atom jumps from excited to ground 
state, depending on the sign of s. 

We will here briefly explain what we mean by a large deviations principle 
and introduce some of the underlying theory. We refer the reader to |16) for 
a complete reference; see fi7| for a comprehensive overview and |18J for an 
introduction to large deviations in the context of statistical mechanics. Large 
deviations is a framework for studying certain types of behaviour of random 
processes (that is, sequences of random variables). In particular, it serves to 
quantify what very rare events are, and to analyse the probabilities of such 
events. Asymptotically (in the sense of the sequence of random variables), such 
rare events will have a probability which decays exponentially, with the rate of 
exponential decay given by a rate function. One of the main results, which we 
will use, is the Gartner-Ellis theorem, which relates this rate function to the lim- 
iting behaviour of the moment generating functions associated to the random 
variables. 

Since a sequence of random variables is entirely characterised by their re- 
spective probability distributions (or laws), the following results are all formu- 
lated in terms of a sequence of probability measures. In particular, the large 
deviation frincvple (LDP) will characterise the limiting behaviour, as n — > oo, of 
a family of probability measures (|^ n )neN on some topological space (X, 23) in 
terms of a rate function. A function I : X — » [0, oo] is called a rate function if it 
is lower semicontinuous (that is, its level sets (x e 1 : I(x) ^ a) are closed); if in 
addition its level sets are compact, we call I a good rate function. The domain of 
I is the set of points in X for which I is finite. 
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The informal idea behind the large deviations principle is the exponential 
decay of the probabilities assigned by the family ( |x n )neN to events xel, which 
we may illustrate by writing 

Mdx) « e- nI(x) dx 

where I is the rate function associated to the family of probability measures. We 
will now state the large deviations principle rigorously. The limiting behaviour 
of the probability measures {\x n } is characterised in terms of asymptotic upper 
and lower bounds on the values that \i n assigns to measurable subsets Pel?. 
The sequence of probability measures {|J. n } satisfies a large deviation principle 
with a rate function I if, for all P € 23, 

- inf I(x) ^ liminf i log (^(P) < limsup ^ log (J. n (P) < — inf I(x). 

xer° n-voo ^o,, xer 

The sentence "|Xn satisfies the LDP" is used as shorthand for satisfies the 
large deviation principle with rate function I." 

We want to prove the LDP for the counting process A t of the atom maser; 
we will do this not by showing that A t satisfies the above definition directly, 
but by applying the Gartner-Ellis theorem to A t . This theorem gives sufficient 
conditions on the sequence of probability measures in order to satisfy an LDP; 
its utility lies in its generality, as it does not require the probability measures to 
be independent. 

Theorem 1 (Gartner-Ellis theorem [16|, pp. 44-55). Let (Z n ) be a sequence of 
random variables in R d with laws \x n . Suppose that the logarithmic moment 
generating function 

A e R d 

exists as an extended real number and is finite in a neighbourhood of the origin, 



A(A) = lim ^logE 

n— foo 
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and let A* denote the Fenchel-Legendre transform of A, given by 

A*(x) = sup {(A. x) — A(A)} . 

AeK d 

If A is an essentially smooth, lower semicontinuous function (e.g. A is differen- 
tiable onl d ) then (Z n ) satisfies a LDP with good rate function A*. 

The large deviation approach involves the use of the perturbed generator 
introduced in Eq. flO) , and its connection to the moment generating function of 
the process A t given by Eq. (IT) . In this setting, the Gartner-Ellis theorem tells 
us that the sequence of atom counts A t satisfies a large deviation principle, 
meaning that 

lim ±logP(A t s$ a) = 1(a). 

t— >-oo 

where 1(a) is the Legendre transform of the function A(s) defined by 

Ms) = Urn l - log E (e sAt ) , (14) 

whenever this function exists and is differentiable. Therefore, to prove the de- 
sired large deviation result we show that this function a>o(s) indeed satisfies 
these conditions. We will do this by showing that o>o(s) is exactly the spectral 
bound of the generator £ s , that is, the largest real part of the eigenvalues of 
£ s . An intermediate step is relating a>o(s) to the spectral radius of the operat- 
ors exp(t£ s ) in the semigroup generated by £ s , which involves using a Krein- 
Rutman theorem to show that this spectral radius is an eigenvalue of exp(t£ s ). 
We will also show that u>o(s) is not only differentiable as a function of s, but 
even analytic. This serves to prove that there always exists a non-zero spectral 
gap between cu s and the rest of the spectrum of £ s ; this spectral analysis is il- 
lustrated in Figj2] Before we go on to state and prove our main result, we need 
to introduce a few more concepts from operator algebra, including the afore- 
mentioned Krein-Rutman theorem. 

Spectral properties of operators 

One of the main results we rely on in proving this large deviations result is 
the Krein-Rutman theorem. In the context of operators on a Banach space, this 
says that the spectral radius of an operator is an eigenvalue with a positive ei- 
genvector, so long as the operator satisfies certain properties. Before we state 
the Krein-Rutman theorem, we introduce some concepts from operator algebra; 
for details we refer the reader to ] 19 ] . 

We are considering a semigroup of maps on the (infinite-dimensional) Banach 
space 23 ( f) ) generated by £ s defined in Eq. (TO) . In our proof we rely on spectral 
properties of such maps. Recall that the spectral radius r(T) of an operator T on 
a Banach space X is given by 

T (T) =sup{|A|:Ae <r(T}} 
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where cr(T) denotes the spectrum of T. The operator T is called compact if the im- 
age of any bounded subset of X is a relatively compact subset; however, we will 
not use this definition directly and will rely on sufficient conditions to obtain 
compactness of the semigroup. The spectrum of a compact operator is always 
a countable set containing 0, that is, if T is compact then 

ff(T) = {0} U {A n : n e N} 

where limn-^ = 0. The resolvent operator R(z, T] is an important tool used to 
relate the spectral properties of the generator L s to that of an element in the 
semigroup; it is defined by 

R(z,T) = (T-p) _1 for all ze C\cx(T). 

In particular, it is used to establish a spectral mapping theorem (see |20|) which 
gives the spectrum of exp(t£ s ) in terms of the spectrum of & sr 

(T(e tL *) = {e tA :AecT(L s )}. 

For this we rely on compactness of the resolvent R(z,L s ), which is a property we 
will prove below. The Krein-Rutman theorem gives more precise information 
about the spectral radius of a compact operator. 

We want to apply a relevant version of the Krein-Rutman theorem to the 
semigroup generated by £ s . Generally such statements are concerned with op- 
erators which preserve the positive cone in a Banach algebra; in particular, to 
obtain uniqueness of the spectral radius as being the only eigenvector of T s (t) 
we need the semigroup to be strongly positive (or positivity improving). This 
means that, given any positive element x e L 2 (f)), we have that T s (t)(x) is 
strictly positive (that is, lies within the interior of the cone of positive elements). 

Theorem 2 (Krein-Rutman |2l) , Ch. 1). Let X be a Banach space, K c X a soid 
cone and T : X — > X a compact linear operator which is strongly positive, that 
is, Tx G int(K) if x G K \ {0}. Then 

1. The spectral radius r(T) is positive and r(T) is a simple eigenvalue with 
an eigenvector x G int(K); there is no other eigenvalue with a positive 
eigenvector; 

2. For any other eigenvalue A ^ r(T) we have |A| < r(T). 

We want to use this theorem to establish that for each t > 0, the spectral 
radius of T s (t) is a simple eigenvalue and T s (t) has no eigenvalues with the 
same modulus. 

3 Results and Methodology 

Our main result can be summarised as follows. 

Theorem 3. The counting process A t satisfies a large deviations principle. 
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This means that, although the behaviour of the spectral bound A(s) is cer- 
tainly interesting (as we will discuss below), there are no non-analyticities; there- 
fore, there are no quantum dynamical phase transitions in this sense. 

The main result is based on showing a Gartner-Ellis theorem for the count- 
ing process A t . This consists of proving a Perron-Frobenius type result for the 
transition operators associated to the modified generator £ s , that is, that the 
spectral radius r s (t) is an eigenvalue of the operator e t£ s , and no other eigen- 
value has the same modulus. The remainder consists of showing that the rate 
function r s (t) is an analytic function of s. The analysis of the spectral properties 
of the operator L s involves choosing a suitable representation for the action of 
the generator, which allows us to show that the semigroup it generates is im- 
mediately compact. 

The steps in the proof are outlined as follows: 

1. Define the L 2 -embedding, following |10 [; the semigroup (T s (t))t^o gives 
rise to a semigroup (T s (t)) t >o on the space of Hilbert-Schmidt operators. 
This new semigroup is used to compute the desired expectation values, 
and we will focus on its spectral properties. 

2. Restrict the associated semigroup (T s (t)) t >o to the commutative subal- 
gebra (of multiplication operators) generated by the number operator N ; 
this defines the generator L s acting on l°°(N). 

3. Show that, in this representation L is self-adjoint and the semigroup gen- 
erated by L is immediately compact. 

4. Use bounded perturbation theory to show that the semigroup generated 
by L s is compact. 

5. Show that the semigroup generated by L s is strictly positive. 

6. Krein-Rutman theorem: the spectral radius r s (t) of T s (t) is an eigenvalue 
with strictly positive eigenvector. 

7. Analyticity of perturbation implies A(s) is analytic. 

Since the transition operator T s (t) is compact for all t > 0, it has discrete 
spectrum; in particular, there is a spectral gap between the top eigenvalue r s (t) 
and the next eigenvalue. Along with the strict positivity of the eigenvector asso- 
ciated to r s (t), we can prove the desired convergence of the moment generating 
functions; analyticity (in particular differentiability) of the rate function com- 
pletes the large deviations principle. 

The first step in the proof is to define a semigroup (T s (t)) t ^ associated 
to (T s (t)) t ^o; this follows fi0| where this is used to find a representation of 
quantum dynamical semigroups. We emphasise that we do not strictly consider 
this new semigroup as a representation of (T s (t))t>0/ rather, it is a separate en- 
tity (T s (t)) t ^o whose expectation values and asymptotic properties coincide. 
We define the new semigroup on the space L 2 ( f) ) of Hilbert-Schmidt operators 
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on the Hilbert space f), which is particularly useful since on this space the gener- 
ator Lo is self-adjoint. This means that we can consider L s as a perturbation of a 
self-adjoint operator, which makes the subsequent analysis significantly easier. 
It is in this setting that we will work for the remainder; we simply refer to L s for 
the associated generator. 

This semigroup (T s (t) ) t ^o leaves the commutative subalgebra generated by 
the number operator N invariant. In the remainder of this article we will con- 
sider the restriction of the semigroup to this diagonal algebra; by T s (t) and L s we 
will implicitly mean the restrictions to this subalgebra. 

As our aim is to apply the Krein-Rutman theorem to the operators T s (t), 
we show that these operators are compact and strictly positive. We show that 
the semigroup generated by L s is an immediately compact semigroup. In or- 
der to show this, we first establish that the generator Lo has compact resolvent 
and generates a semigroup of compact operators. Then, upon showing that the 
perturbation term in L s = L + A s is bounded, we apply bounded perturbation 
theory to obtain the desired result. 

Because T s (t) is compact, we can apply the Krein-Rutman theorem for com- 
pact operators to claim that the spectral radius r s (t) — exp(tA(s)) is an eigen- 
value of T s (t), which is associated to a positive eigenvector |R S ). Furthermore, 
because T s (t] is positivity improving, the eigenvector |R S ) is strictly positive. 
This then allows us to conclude that the convergence required for the Gartner- 
Ellis theorem holds, that is, 

lim -logE(e sAf ) = A(s) 

t->oo t 

which establishes the large deviation principle. Finally, because L s is an analytic 
perturbation of L, the function A(s) is itself analytic. 

4 Details of proof 

The proof of our main result starts with a study of the spectral properties of the 
perturbed generator L s . Before we begin, we define an associated semigroup 
of the dynamics by embedding the algebra 23(f)) in the Hilbert space L 2 (f)) of 
Hilbert-Schmidt operators (we adapt this technique from JT0)). We will refer to 
this associated semigroup as an embedding of the original semigroup. 

Recall that X £ 23(f)) is called a Hilbert-Schmidt operator if it has finite 
Hilbert-Schmidt norm, that is, 

||X||2j„ s =Tr|X| 2 < oo. 

The space of all Hilbert-Schmidt operators, denoted L 2 (f)), is a Hilbert space 
with inner product 

(X, Y) H _ S = Tr (X*Y) for all X. Y e L 2 (f)). (15) 

This map which embeds 23(f)) into L 2 (f)) is defined in terms of the stationary 
state p ss (see Eq. fl2) ) of the generator L. 
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Following [TO) , we define the embedding i as 

i:23(f)WL 2 (f)) 
: x i ^ PsWsY 4 - 

Using this map i, we can define the embedding of operators on 23 ( t) ); with an 
operator T : D{7) 23(f)) we associate an operator T on L 2 ([)) by 

T(i(X)) = i(T(X)) for all X e T>[7) c 23(f)). 

This embedding is used to define the associated semigroup; we formally define 
the generator L s on L 2 (f)) by 

L s (X)=i(£ s (i- 1 (X))). (16) 
It will be helpful to write L s in the perturbative form 

L S (X)=L(X)+A S (X) (17) 



where, recalling Eq. 10 L is the generator obtained by applying the embedding 
to the original generator L and A s is the perturbation arising from biasing the 
generator. By restricting to finite-rank elements of L 2 ( f) ) it is impossible to find 
an explicit expression for the action of L s . 

We will be using this associated semigroup (T s (t)) t ^o m the remainder of 
this paper. To prove the large deviations principle for the counting process A t , 
we express the expectation value in Eq. (TT) in terms of an expectation value on 
L 2 (f)). Using Eq. (16) we are able to show that 



Tr (p in T s (t) (1) ) = Tr ( Pln T s (t) ( p, 1 / 2 )) 



where we have defined p in = P ss 174 pinPss 174 ; in terms of the inner product on 
L 2 (f)) defined in Eq. (15) , this expectation value reads 



Tr(p in T s (t)(l)) = (p in ,T s (t)( Ps 1 / 2 ) / , 

\ / H-S 

Here the action of the semigroup generated by L s is still taken in the Heisenberg 
picture. 

The semigroup (T s (t)) t ^o leaves the commutative subalgebra of multiplica- 
tion operators invariant; in particulari the stationary state p ss is a multiplication 
operator. In the remainder of this paper we will only be considering the action 
of L s on the commutative subalgebra N of L 2 ( f) ) consisting of multiplication 
operators. 

In particular, for any f G D(L S ) we have 

L s (M f ) = M q 
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where g € l 2 (N) is given by 

g(n] = A n _^ n f(n - 1) - (A 2 + ^)f(n) + A n ^ n+1 f(n + 1) (18) 
+ (e- 1 -l)^^N ex sin 2 (4)Vn~TT)f(n+l). (19) 

A n 

Here the rates A^, u.^ are given by 



A k = ^/v(k+l) + N ex sin 2 (cf>Vk+T), 

In the unperturbed case s = 0, Eq. (18) shows that the restriction of L s to the 
commutative subalgebra of L 2 ( f) ) generated by I 2 (N) corresponds to a birth-death 
process from classical probability theory with rates A^, u^. Therefore L coincides 
with the generator of this classical birth-death process, where for any f £ D(L) 

L(f)(n) = A n _innf(n- 1) - (A 2 , + ^)f(n) + A n ^ n+1 f(n + 1). (20) 

The perturbation in Eq. (IT) takes the following explicit form: for any f g D[L) 

A s (f)(n) = (e- 1 -l)^±iN ex sin 2 (cl 5 Vn+l)f(n+l). (21) 

A n 

This defines the representation L s on the Hilbert space l 2 (N) of the perturbed 
associated generator; abusing notation, we clearly have 

L s (M f ) = M Ls(f) for all f e D(L). 

The remainder of the proof will focus on the spectral properties of the generator 
L s . 

Lemma 1. The semigroup generated by L is immediately compact. 

Proof. Since the L 2 -embedding defines L as a self -adjoint operator on the Hilbert 
space l 2 (N), the spectral theorem applies. Let the spectral decomposition of L 
be of the form 

L = Y_ A nPn, 

n^O 

where cr(L) = (A n ) n ^o is the spectrum of L and P n is the eigenprojection asso- 
ciated to A n . Then the spectral theorem for self -adjoint operators states that 

T t = Y_ e tA -P n for all t > 0. (22) 

n^O 

This ansatz will be used to show that T t is compact for every t > 0; in particu- 
lar, we will show that T t may be approximated in norm by the partial sums in 
Eq. p2) . In order for this to be valid, we show that L has discrete spectrum (in 
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the sense of having a countable number of eigenvalues and finite-dimensional 
eigenspaces), and the eigenvalues of L diverge to — oo. 

In order to establish that L has discrete spectrum, we follow the argument 
used in 1 22]. Note that the action (see Eq. ( (20) ) of L on l 2 (N) maybe represented 
by the infinite-dimensional tri-diagonal matrix 



L = 



f -(A 2 , + Aom \_ 

o Am 2 -m + v%) 



(23) 



acting on infinite column vectors representing f G l 2 (N). We split L up into a 
sum 

L = A + B 

of two infinite-dimensional matrices; A containing the diagonal elements in 
Eq. |23|, and B containing the off-diagonal elements. Then L can be regarded 
as a perturbation of A; in particular, the perturbation is relatively bounded in 
the sense of Kato |23 J, that is, 

j|Bx|| sC ||Ax|| for all x G D[A). 

The resolvent of A is compact, and using Kato's result on relatively bounded 
perturbations (p. 214) we find that the resolvent of L = A + B is also compact. 
Since the self -adjoint operator L is closed (see |24|, p. 256), it follows (Kato, p. 
187) that the spectrum of L is discrete. 

Consider the the partial sums of the exponential of the generator L defined 

by 



T< N '(t) = X 



n=0 

where {Ai, A 2 , . . .} is the spectrum of L and P n is the projection onto the eigen- 
space corresponding to A n . The real eigenvalues {Ai, A 2 , . . .} are bounded from 
above and by using the spectral mapping theorem on the spectrum of the com- 
pact resolvent of L we conclude that A n — > — oo. Therefore e tL may be approx- 
imated by the finite-rank partial sums and e tL is compact for all t > 0. □ 

Having established that the semigroup generated by L is compact, we use 
a bounded perturbation argument to extend this property to the semigroup 
generated by L s . 

Corollary 1. The semigroup generated by L s is immediately compact. 

Proof. Recall that the perturbation in Eq. < [T7) may be represented on the com- 
mutative subalgebra l 2 (N) as 

L s = L + A s 

where A s is defined in Eq. | [2~T| . Since A s is a bounded operator, the semigroup 
generated by L s is also immediately compact (see 1 13], p. 167). □ 
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The Krein-Rutman theorem also requires that the operators T s (t) are all 
strictly positive, that is, (en,|T s (t)(p)|e n ) > for all n ^ 0. This follows from 
a Dyson expansion (see [25 1) applied to the generator L s = £° + L s , where L° 
is the generator obtained from L s by taking 

£ s (p)=£ (p) + e- s L lP q + L 2 pq. 

The resulting generator L° is associated to a classical birth-death process and 
so the commutative restriction is irreducible, and so strictly positive. As the 
leading term in the Dyson expansion involves only L°, and all other terms (in- 
volving the jump operators I_i , L2) are positive, this suffices to conclude that the 
operators T s (t) are strictly positive. 

Our goal is to show that the spectral bound A(s), defined as the largest real 
part of the spectrum 

A(s) ={Re(A) :Ae <r(L s )} 

plays the role of the large deviations rate function as introduced in Eqs. (TT),([T4). 
The first step in establishing this follows from compactness of the operators 
T s (t) and the Krein-Rutman theorem. 

Corollary 2. The spectral radius r s (t) of T s (t) is an eigenvalue of T s (t) and has 
associated to it a strictly positive eigenvector. 

Proof. Compactness of the semigroup (T s (t)) t ^ leads( [13 1, p. 281) to a spectral 
mapping theorem of the form 

e ttr(L s ) = ff (T s (t))\{0} for all 0. 

Since each operator T s (t) is compact and strictly positive, we can apply the 
Krein-Rutman theorem for compact operators (see Thm. |2]above). □ 

Since r s (t) e ff(T s (t)), using the spectral mapping theorem this means that 

r s (t) = e tA(s) 

where A(s) is the spectral bound of L s ; any other eigenvalue A £ cr(T s (t)) satis- 
fies |A| < r s (t). 

To complete the proof of the large deviations principle, we need to show that 
the spectral radius r s (t) of T s (t) is a differentiable function of s. Using analytic 
perturbation theory for the generator L s , we can in fact prove that the spectral 
bound s(L s ) is an analytic function of s. 

The results used on analytic perturbation theory are found in ] 26 1 . Given a 
closed operator T' we define the family of perturbations T(s) by 

T(s) =T'-sH, 

where H is a closed operator with the same domain as T'. The main result 
(Prop. 3.25, p. 141) on analycity is that any isolated eigenvalue (of finite mul- 
tiplicity) and its associated eigenprojection are analytic functions of s in some 
disc around s = 0. We apply this result to the family of perturbations 

L s =L + A S , 
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Figure 3: Comparison of classical birth (solid) and death (dashed) rates (see 
text); stationary distribution 

for which it is clear that the spectral bound of L s is an analytic function of s. 
Since, by the Krein-Rutman theorem, the spectral bound of L s and the spectral 
radius r s (t) are related by 

r s (t)=e s < L *\ 

we immediately see that r s (t) is also an analytic function of s. As a consequence 
of this, the spectral radius remains isolated when viewed as a function of s 
(which should not be confused with its isolation for any fixed s since that trivi- 
ally follows from compactness of the semigroup). 

5 Numerics and discussion 

We have shown above that the counting process A t satisfies a large deviations 
principle; in particular, the spectral radius r s (t) is an analytic function of the 
large deviations parameter s. We will now briefly consider the interpretation 
of the large deviations rate function in the context of quantum dynamical phase 
transitions. 

In [1 J the interpretation of the large deviations rate function as an indicator 
of quantum dynamical phase transitions. Non-analyticities in the rate function, 
such as discontinuities in the first derivative correspond to the spectral gap be- 
coming very small (see Fig.[2j. At these values for <x we also see that the station- 
ary state becomes bistable; that is, instead of a single peak in the mean photon 
number distribution, we see two separated peaks (see Fig. [TJ. This means that 
the state of the cavity can be in either of two bistable phases, one more active 
than the other; these two phases become apparent in simulations. The num- 
ber of excitations in the cavity spends some time in a more active region before 
making a transition to a less active region. The occurrence of this behaviour for 
some values of a signifies a different dynamical region from other values of oc, 
and this is a quantum dynamical phase transition. 

Now consider the jumps in the mean photon number (solid line in Fig. |3|. 
The values for oc associated to the midpoint of the jump (e.g. <x c « 6.7) are 
associated to points where the derivative A'(s) of the spectral bound A crosses 
the s-axis in the following sense. The derivative A'(s) shows clearly different 
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l"(s) for counling measurement with N ^ = 100 




(b) Crossing at <x c for N ex =(c) Mean photon num- 
(a) A'(s) (asinFig.|2J 20, 30, 90, 100 (varying s axes) ber (see text) 



Figure 4: Analysis of crossing near a c 

values when comparing either sides of the 'phase lines' visible in Fig. [4] (a). 
Although this figure suggests that this phase line approaches the s = axis 
strictly from s < 0, Fig. |4] (b) shows that this line crosses over to the s > 
region. The angle at which the crossing occurs increases with N ex (note the 
changing s ranges: the top two figures are with |s| ^ 0.01 while the bottom two 
figures are with |s| ^ 10~ 8 ). 

We are interested in the s = behaviour since, if A'(0) is discontinuous, this 
means there is a dynamical phase transition: the mean (A t ) (first moments) 
depends on whether the derivative is taken from the left or right. Since the 
phase line does cross the s = axis, this does not necessarily preclude such a 
phase transition. However, we know mathematically that the spectral bound A 
is an analytic function of s, so a phase transition cannot occur. 

Accuracy of the numerical derivative of the rate function is illustrated in 
Fig. [I] (c). If the derivative with respect to s is taken with the differential As 
small enough (10~ 8 ), the derivative follows the mean photon number exactly. 
For larger As, the derivative follows the 'wrong' phase and thus results in the 
wrong mean being computed, which is clearly visible in the top figure. 

Another way of looking at these metastable phases is by considering the 
classical birth-death process associated to the quantum dynamical system. Fol- 
lowing the discussion in p7| , we consider the values of a for which the 'birth 
rates A^ exceed the 'death rates [i^ along with the dynamical nature of these 
crossover points. In Fig. [3] we see these crossover points; for example, at a = 1 
the rates cross at a maximum for the death rate, which in the stationary state 
corresponds to the sharpest increase in the mean photon number. At a = 4.6 
the rates intersect in more than one location for the first time, which is reflec- 
ted in the photon number distribution as the start of the bistable region. (Note 
that this figure is independent of the variable N ex , as x here is an appropriate 
rescaling of n in the limit N ex — > oo.) 
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6 Concluding remarks 



We have studied the counting process associated to photon number trajectories 
of the cavity state in the atom maser, and shown that this process satisfies a 
large deviations principle. Our approach could be extended to measurements 
other than the counting of ground-state atoms used here, for example, homo- 
dyne detection. This means taking a different perturbation of the generator to 
analyse the large deviations behaviour. However, the resulting dynamics will 
no longer have the nice property of leaving the diagonal subalgebra invariant, 
so the algebra and the numerical analysis will be more involved, requiring non- 
commutative analysis. The method could also be extended to any immediately 
compact, strictly positive semigroup; that is, the exact form of the generator 
does not directly affect the principles used here. Related results on a general- 
isation large deviations behaviour of the counting statistics in quantum spin 
chains (as addressed in |7|) will be discussed in a forthcoming publication. 
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A Equality of expectation values for L s and L s 

Instead of considering the embedded generator L s as a representation of L s on 
L 2 (f)), we should interpret these as generators of separate semigroups, one on 
L 2 (rj) and the other on 23(f)), and show that they agree in the expectation value 
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we are interested in. That is, show that the quantities defined by the moment 
generating function 

E(e sAt ) =Tr( Pin T s (t)(l)) 

and the inner product 

Tr(p in e tL *(p 1/2 )) = (Pin,Ts(t)(p 1/2 )) H -s 

are finite, and equal, for all finite-rank initial states and all t > 0. Note 
that the moment generating function should be finite since it is bounded by the 
moment generating function for a Poisson process 



E(e sAt ) =^e sk P p (A t = k) (24) 

k^O 

< Y_ e sk Pp oiS son(k) (25) 



k>0 



since the arrival times of the atoms are governed by a Poisson distribution. 

Following a similar argument in |28], we show that there exists a time T > 
such that for all t < T, the series expansion 



k! 

k^O 



defines a bounded operator for all finite rank operators x e 23(f)). In the lan- 
guage of operator analysis, all finite-rank operators are analytic vectors for L s . 
This is done by considering rank one operators of the form |e n ) (e m | and finding 
a T > such that the series 



V- 



\£>s[\e n )(er 



k! 

k>0 

converges for all t < T. Such a T is readily found (see |28|) and depends on the 
coefficients \x, A appearing in the master equation. 

Having shown that all finite rank operators are analytic vectors for L s , we 
continue by noting that given any operator x € S ( f) ] we can define its (finite 
rank) projection 

x n = P n xP n 

where P n projects onto the subspace of rank n operators. We are interested in 
the action of the semigroup on trace-class operators; the generator of this p 're- 
dual semigroup is denoted by £ s>t and is induced by L s by the duality relation 

Tr(e t£s *(p)x) =Tr(pe t£s (x)) 

for all p € Ti([)) and all x € T>{t)). 
We want to show equality in 



Tr(p in e t£ *(l)) =Tr(p in e tL *(p 1/2 )) 
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for all t > 0, (for now) for an appropriately chosen set of initial states Pi n . Par- 
tition the time t by writing t = ti + . . . + t\ where < T. The main idea 
behind showing the above equality is by applying sequentially the analyticity 
of finite-rank operators and weak-* continuity to 'move over' each element of 
the semigroup e ti£ to the trace-class operators and then the Hilbert-Schmidt 
operators, and repeat this until the entire semigroup lives on the other side. 
Using the projections P n introduced above, we find that by weak-* continuity 
and the semigroup property we have 

Tr (pe (tl+ - +tl)£s (l)) =^1^ Tr (pe tl£s (P n e (t2+ - +tl)£s (l)P n )) . 

Define the finite-rank operator A n = P n e (t2+ - +tl)£s (l)P n ; then since this 
is an analytic vector by the above argument, we may express the RHS as a power 
series and thus we obtain 

\ k>0 ' / \ k>0 / 

by definition of the embedded generator L s . Write the embedded initial state 

as p = p7s 1/4 pp7s 1/4 and this limit may be written as an inner product in the 
Hilbert-Schmidt space as 

lim Tr (pe tlL * (pi/ 4 P n AP n pi/ 4 )) = lim (p, e tlL * (A n )) (26) 



lim (e tlL »(p),A n ). (27) 

n— too 



It now suffices to show that A n -> A weakly; this is clear, since for any x e 7i [ f) ) 
we have 

(t, P n p^ 4 Api/ 4 P n ) = Tr (TP n pi/ 4 Ap 4 / 4 P n ) (28) 
= Tr (p 4 / 4 Tpi/ 4 P n AP n ) (29) 
Tr (p^ 4 xpi/ 4 A) ; (30) 



whenever Ps{ 4 Tpss 4 € Ti(f)), 



Tr(pV 4 Tp 4 / 4 )=Tr(p 4 / 2 T)=(p 4 / 2 ,T)<oo 



so 



P n pj{ 4 Apj{ 4 P n — > p 4 { 4 Apg{ 4 
weakly. Therefore the limit becomes 

Tr(pe^ L '(pi/ 4 e^ + - +t ») £ -(l)pi/ 4 )) (31) 
= lim Tr (e^ L » (p^V 2 ^ (P n BP n ) p 4 { 4 ) . (32) 
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Similar to before, we define B n = P n BP n where B = e (t3+ - +tl)£s (1). By re- 
peating the above argument until all the t t have been exhausted we arrive at 
the desired expression 

Tr (e^ + ■ +tl)L » (p)pi/ 2 ) = Tr (pe tL * (p^ 2 )) . 

This shows that the semigroup L s on the Hilbert-Schmidt space gives rise to the 
desired expectation values. 
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